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[1] Numerical Analysis via DERIVE, Steven Schonefeld, 530 pages 
MathWare, 604 E. Mumford Dr. , Urbana, IL 6 1 80 1 , 1 994 

[2] Precalculus Investigations using DERIVE, David M. Mathews, 128 pages 
Harper Collins College Publishers, 1994 

[3] The DERIVE Calculus Workbook, L. Townsley Kulich and Barbara Victor 
Harper Collins College Publishers, 1994 

[4] Mathematics with Excel, David Sjostrand, 189 pages 
Chartwell Bratt, 1994 

[5] Matematikk for 0konom og Samfunnsfag, H.Bjomestad and others, 489 pages 
(Mathematics for economics and social studies; the first textbook of mathematics using DERIVE in 
Norway)Hoyskole Forlaget - Norwegian Academic Press, 1994 

[6] LAB Projects, NSF Calculus Institute, Using Computer Algebra Systems, 510 pages 
Editor: J.Birdsall, Santa Fe Community College, Gainesville, FL 32606 

[7] Derive per la scuola, S. Rossetto, 138 pages 
McGraw - Hill Libri, Milano/Italy, 1992 

(Thank you Giuseppe, the other book you mentioned is already on our shelf!) 


e-mail .... e-mail .... e-mail .... e-mail .... e-mail .... e-mail .... e-mail 


We have received many requests concerning e-mail numbers of Soft Warehouse or of the BBS 

I am glad to be able to give an answer now. Yes, there are some e-mail numbers: 

Soft Warehose: 74521.3051 (CompuServe) 

74521 .3051@compuserve.com (Internet) 

DERIVE-NEWS: Derive-News is a new electronic discussion list set up after the International De- 
rive Conference at Plymouth in July. You can join by sending a message to 

mailbase@mailbase.ac.uk with the single line 

join derive-news firstname lastname 

(eg join derive-news pam bishop) as the text of the message. 

European DERIVE Bulletin Board: cain@can.nl (Internet) 

e-mail .... e-mail .... e-mail .... e-mail .... e-mail .... e-mail .... e-mail 


For our German members: 

In Deutschland sind viele Schulen an das „Offene Deutsche Schulnetz" ODS angeschlossen, das den 
Schulen kostenlosen Internet-Zugang ermoglicht. Dort gibt es sogenannte NEWS-Gruppen fur gewis- 
se Themen - so zB die Gruppe “schule.math". (Hinweis von G.Noll, Erpel/Rhein) 

H. Scheuermann wurde gerne mit anderen Kollegen in Kontakt treten und bittet mich seine Adresse 
zu veroffentlichen. (Aus Griinden des Datenschutzes will ich keine Liste der Mitglieder weiterleiten.) 
Also, wer Lust auf einen Gedankenaustausch hat wendet sich an: Hellmut Scheuermann, Ubierstr. 1 3, 
65719 Hofheim/Taunus, Tel.: 06192/27566. 


Two offers: 

I have a copy of DERIVE, Version 3 Release Notes and I have a copy of “The Evolution of 
DERIVE“. If you are interested in one of these papers (- or both) then please write or call I'll 
send you the papers. 
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Liebe DUG-Mitglieder, 


Dear DUG Members, 


Ein schoner Sommer - fur Osterreich und einen 
groBen Teil Europas - geht zu Ende und ich freue 
mich, Ihnen die 15. Ausgabe des DNL prasentieren 
zu konnen. Fur die DERIVE-Gemeinde war die 
Konferenz in Plymouth sicher das herausragende 
Ereignis dieses Sommers. Hier ist zu wenig Platz, 
jeden einzelnen der vielen hervorragenden Vortra- 
ge zu wiirdigen. Sie konnen sich aber sicher vor- 
stellen, daB besonders der Hauptvortrag von David 
Stoutemyer und Albert Rich iiber das „Innenleben“ 
von DERIVE und iiber die Moglichkeiten von 
DERIVE 3.x ganz besonderes Interesse gefunden 
haben. An dieser Stelle mochte ich im Namen der 
DUG John Berry und seinem Team fiir ihren Ein- 
satz recht herzlich danken. 

Ich habe im letzten DNL einen Beitrag von Sergey 
Biryukow versprochen, der die Beschriftung und 
Skalierung von DERIVE-Plots aus der DERIVE- 
Umgebung heraus ermoglicht. Da ich Ihnen nicht 
zumuten will, in muhsamer Arbeit die Daten fiir die 
Zeichensatze einzugeben, werde ich diesen Artikel 
in der nachsten Nummer drucken und kann die 
zugehorigen Datendateien auf der Diskette 94 
gleich mitliefem. Bis dahin bitte ich um Geduld. 
(Wenn jemand wirklich nicht mehr warten kann, 
dann soil er mich bitte anrufen.) Sergej iiberraschte 
uns in Plymouth nicht nur mit einem spannenden 
Vortrag iiber DERIVE-Anwendungen (Optik) 
sondern auch mit einem DERIVE - Gedicht (Seite 
48). 

Auch die Freunde von Th. Weths „Kurvenlexikon“ 
bitte ich um Geduld fiir den nachsten DNL. Aus 
Platzgriinden einerseits und aus Riicksicht auf un- 
sere vielen nicht deutschsprachigen Freunde ande- 
rerseits wird die Serie mit Folge 5 (mit der Kon- 
choide des Nikomedes) im nachsten DNL fortge- 
setzt. 

Sie fmden dieses mal auf unserem Biicherbrett auch 
ein Buch, das sich nur am Rande mit DERIVE 
beschaftigt. David Sjostrand ist ein sehr aktives 
Mitglied der DERIVE-Gruppe und schrieb ein 
Buch iiber Mathematik mit Excel. Fiir ihn - und fiir 
viele andere - sind Tabellenkalkulation und 
DERIVE kein entweder - oder sondern ein sowohl - 
als auch. 

Aus dem angekiindigten Splines - Wettkampf 
Bohm - Reichel hat sich durch einen Beitrag von 
LeoKlingen, Bonn, ein Dreikampf entwickelt. Nach 
meiner ersten Einschatzung diirfte die reiche Erfah- 
rung Dr. Klingens ihm einen betrachtlichen Vor- 
sprung sichern. 


Mit den besten GriiBen 
bis zum nachsten Mal 



A wonderful summer - for Austria and a large part 
of Europe - is going to have an end now and I have 
the pleasure to present the 15th issue of the DNL. 
For the DERIVE-community the Plymouth Confer- 
ence in July was the outstanding event of this sum- 
mer. Here is not room enough to mention honora- 
bly each one of the interesting lectures. But you 
surely can imagine that we all especially appreci- 
ated David Stoutemyer's and Albert Rich's plenary 
lecture about “DERIVE 's interior 44 and about 
DERIVE 3.x's new features. Let me say a warm 
THANK YOU in the DUG's name to John Berry 
and his team for organizing the conference and 
making it a success. 

In the last DNL Eve announced a contribution from 
Sergey Biryukow which enables us to label and to 
scale a DERIVE-plot within the DERIVE- 
environment. Because it is a very hard and boring 
work to edit the fonts' data I've decided to publish 
this article in the next DNL, so you can find then 
the .MTH-file and the data files on the diskette of 
'94. I hope that is in your sense. Sergey surprised 
us in Plymouth not only with an excellent lecture 
about DERIVE applications (optics) but also with a 
DERIVE-poem 
(page 8). 

I also ask the friends of Thomas Weth's “Lexicon 
of curves 44 for patience. There are two reasons to 
continue his series (with the Konchoide of Niko- 
medes) in the next DNL: arrangement in the layout 
and consideration for our many non German speak- 
ing friends. 

On the Book-Shelf you can find this time one title 
which is dealing with DERIVE only very little. 
D. Sjostrand is a very active DERI VI AN and he 
wrote a book about teaching Mathematics with 
Excel. For him - and for many others - DERIVE 
and spread sheets don't represent a either - or but 
an as well - as. 

The promised spline - competition Bohm vs Rei- 
chel seems to become a triathlon because of a 
SPLINE.MTH file from Leo Klingen, Bonn. I think 
that Dr. Klingen 's rich experience will save him a 
considerable lead. 

Sicerely yours 


ie Tagungsbande von Krems '92 „Teaching Ma- 
thematics with DERIVE 66 und Krems '93 
„DERIVE in Education 66 sind zum Preis von 
oS 330.- /oS 450.- bzw. DM 50.- / DM 65.- (inkl. 
Verpackung und Porto) bei mir zu beziehen. Anruf, 
Fax oder kurzes Schreiben geniigen. 
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The DERIVE-NEWSLETTER is the Bulle- 
tin of the DERIVE User Group. It is pub- 
lished at least four times a year with a 
contents of 30 pages minimum. The goals 
of the D-N-L are to enable the exchange 
of experiences made with DERIVE as well 
as to create a group to discuss the possi- 
bilities of new methodical and didactical 
manners in teaching mathematics. 


Preview: (Contributions for the next issues): 

Stability of systems of ODEs, Kozubik, SLO 
Prime Iterating Number Generators, Wild, UK 

Graphic Integration, Probability Theory, Linear Programming, Bohm, AUS 

DERIVE in Austrian Schools, some examples, Lechner, Voigt, Eisler a.o., AUS 

Tilgung fremderregter Schwingungen, Klingen, GER 

Der Fermat - Punkt im Dreieck, Geyer, GER 

Continued Fractions and the Bessel Functions, Cordoba a.o., ESP 

Turtle Commands in DERIVE, Lechner, AUS 

DREIECK. MTH, Wadsack, AUS 

Newton Method and Ill-Conditioned Problems, Lopes, POR 

2D Plots Labeling, Biryukov, RUS 

Plotting t-periodic functions, Verhoosel, NED 

IMP Logo and Misguided Missiles, Sawada, HAWAII 

“Reverse" Discussion of Curves, Reichel, AUS 

Reichel - Klingen - Bohm - Splines, A thriathlon, AUS & GER 

3D Geometry, Reichel, AUS 

Parallel- and Central Projection, Bohm, AUS 

Conic Sections, Fuchs, AUS 

The IMP Logo, Sawada, USA 

Setif, France; Vermeylen, Belgium; Leinbach, USA and others 


Editor: Mag. Josef Bohm 
A-3042 Wiirmla 
D 'Lust 1 
Austria 

Phone: 43-(0)2275/8207 


Contributions: 

Please send all contributions to the Editor. 
Non-English speakers are encouraged to 
write their contributions in English to re- 
inforce the international touch of the 
D-N-L. It must be said, though, that non- 
English articles are very welcome none- 
theless. Your contributions will be edited 
but not assessed. By submitting articles 
the author gives his consent for reprinting 
it in D-N-L. The more contributions you 
will send, the more lively and richer in 
contents the DERIVE Newsletter will be 


Impressum: 

Medieninhaber: DERIVE User Group, A-3042 Wurmla, D'Lust 1, AUSTRIA 
Richtung: Fachzeitschrift 
Herausgeber: Mag.Josef Bohm 
Herstellung: Selbstverlag 
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Steven Schonefeld, Angola, IN, USA 

Can DERIVE decide if two numbers are equal? 

The answer is definitely YES if you are using version 2.58 or later. In some earlier versions, 
DERIVE's IF(,,,) function cannot always decide whether two numbers are equal, especially if one or both 
of the numbers are complex. For example, the DERIVE function COMPARER, y), defined below, should 
decide whether the numbers x and y are "Equal" or "Not equal". 

COMPARE (x, y) := IF (x = y, "Equal", "Not equal", "?") 

COMPARE (3.14, n) [Simplifies to "Not equal" 

COMPARE (1/(1 + i) , (1 - i)/2) [Simplifies to "Equal" 

COMPARE (1, i) [Simplifies to "?" when we use version 2.09. However, 

COMPARE (1, i) [Simplifies to "Not equal" when we use version 2.58. 

There is a simple way to fix this problem so that you get the correct results for all versions of 
DERIVE containing the IF(,,,) function. This fix relies on the fact that numbers x and y are equal if, and 
only if, the absolute value of their difference is zero. Thus, we replace "x = y" with "ABS(x - y) = 0" in the 
function COMPARER) to get. 

COMPARE_2 (x, y) := IF (ABS (x - y) =0, "Equal", "Not equal", "?") 

The function COMPARE_2(x, y) now gives the correct results for all versions of Derive containing the 
IF(,„) function. See DERIVE 6: 


No problems with recent versions of DERIVE 

compared, y) := 

If x = y 

#1: "Equal 11 

"Not Equal" 

"Cannot decide" 

#2: compare(3 . 14 , m) = Not Equal 

f 1 1 - 


#3: compare 


= Equal 


L 1 + i 2 J 
#4: compare(l, O = Not Equal 


Below is an implementation of the secant method (for approximating solutions to F(x) = 0). We 
see that DERIVE needs to determine whether the two numbers ELEMENT(v, 3) (= vl3) and 
ELEMENT(v, 4) (=vJ4) are different in order to prevent division by zero in the calculation of NEW_5(v). 

#1: Ffr) := 

NEW_5(v) := 

If v|3 = v^4 
# 2 : 0 

Vi2-V44/(V43 - vj,4) 

NEW _V(v, 6) := |~v + 6, 6, v , F(v + 5)1 

#3: L 1 4 1 J 


#4: 


SECANT (a, 5, n) := ITERATES (NEW_V(v, NEW_5(v)), v, [a, 5, F(a - 5), F(a)], n) 
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Terence Etchells, Bolton. UK 

Here is a little function for the inverse normal distribution: 




# 1 : 


INVNORMCp, p, 


cr) := ITERATE 


a 


SIGN 0). ERF 


7 2 - a 


2 -cr 


2 



P 


2 2 2 2 2 
- a /(2-cr ) + a-p/cr - |j /(2-cr ) 

^2-e 


a, 


v 


2 ■ Jtt ■ | cr 


Example: 


# 2 : 


#3: 

#4: 


NORMAL (109 . S, 105, 7.5) = 

NORMAL (109. S, 105, 7.5) = 
INVNORM(0 . 73S9 , 105, 7.5) 


ERF 


S-^2 


25 


0 . 73S9137003 
= 109. 7996 S39 


+ 1 


Working now with NSOLVE: 

#5: NSOLVE (NORMAL (sc , 105, 7.5) = 0.73S9, x) = (x = 109.7996S39) 
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Josef Bohm, Wiirmla. A 

Two tipps for working with DERIVE: 

Random Numbers 

I worked in the classroom and wanted to demonstrate an experiment using random numbers. All the 
students loaded the file, simplified and then we compared the results. You will understand that I was 
not very happy when we found out, that all the results were identical. We tried once more: other re- 
sults, but identical again 

So my DERIVE show was not a big success. What to do? 

I took the manual and tried to find a way. That is it: 

Before starting the simulation simplify the command 

RANDOM (0) 

then go on in the usual way. Doing this you initialize the random number generator with the 
time elapsed since starting DERIVE. 


Text processing and DERIVE 

There are problems to convert DERIVE-files saved in extended ASCII-Code, into WORD for Windows 
(in ANSI-Code). Now in WORD for Windows 6.x it is possible: 

• Save the DERIVE - file with Transfer Print File FILENAME.TXT. 

• Open the file from WINWORD and then format with the font MS LineDraw. Compress the lines, 
eg. line distance 0.75. Don't be surprised that the text now may look a bit strange. A result of this 
effort can be seen above. 

(This advice is no more valid in times of DERIVE for Windows.) 
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Message 2622: From HARALD LANG to PUBLIC about BUG IN V 2.52 

Version 2.52 has a bug. Solving a-b*tan(x) A 3 produces a wrong answer! Are there more bugs? 

Message 2625: From DAN to HARALD LANG about #2622 / BUG IN V 2.52 

Could you give your answer? What branch was your default on? Did you try the real branch? 

— Dan Frezza bz223@Cleveland.Freenet.Edu 

Message 2627: From HARALD LANG to DAN about #2625 / BUG IN V 2.52 

With Branch=real, solving e.g. 8-8*TAN(x) A 3=0 produces the answer x=ATAN(l/4). In general, 
a-b*TAN(x) A 3 gives the answer x=ATAN[(a A (l/3)/b] rather than x=ATAN[(a/b) A (l/3)]. Also other 
powers, like 2, produces an incorrect answer. Have You got the same problem? 

Message 2630: From DAN to HARALD LANG about #2627 / BUG IN V 2.52 

Harold, I believe you have found a bug. I tried the same thing, and got the same answer. 

How is DERIVE 6 responding? 


SOLVE (a - b-TAN(x) = 0, x) 


r 

r i/Bfi M'V 


r 1/Bfi ^3. L y 

x = - ATAN 

3 G ■ 2 J 

- 77 Y X = 77 - ATAN 

T 2 ~ 2 J 

1/3 

1/3 




^ b ^ 


^ b ^ 


ATAN 


ATAN 



v x = - ATAN 


r 1/3 r 1 

1- ' 


r 1/3 r 1 

r- ■ 

a 

2 J 

- 77 y x = 77 - ATAN 

a | T - 

2 J 

1/3 

Lb J 

1/3 

Lb J 


v x = ATAN 


1/3 


1/3 


77 v x = ATAN 


1/3 


1/3 


L b 


77 v x = ATAN 


1/3 


1/3 


S0LVE(S - S- TANGO =0, x) = x = 


S ■ 77 OLN(2 - ^3) 


5 ■ 77 L-LNG/3 + 2) 


3 -77 


i. ■ LN(2 - A) 3-77 i ■ LNC/3 +2) 77 l-LN(2 - ^3) 77 l- LN(A + 2) 5-77 

Y X = - Y X = - - Y X = - - Y X = 

2 42 42 42 4 


3 ■ 77 77 ^ 

V X = - Y X = 

4 4 J 

3 

SOLVE (8 - S-TANOO = 0, x, Real) 

5 - 77 3 ■ 77 77 

X = Y X = - Y X = 

4 4 4 


Message 3098: From HARALD LANG to HADUD about #3091 / SIMPLIFYING CUBE 
ROOTS 

Here is another example: 8*SIN(7i/18)*COS(7i/9)*COS(27i/9) 

This is equal to 1 (new challenge: prove it!). We can get reasonably convinced by approximating it by 
DERIVE, but DERIVE fails to simplify it to 1. A Bug!! Hello, Soft Warehouse, are you there <the- 
re...echo...echo> :-) 
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Message 3100: From SOFT WAREHOUSE to HARALD LANG about #3098 / SIMPLIFYING 
TRIG EXPRESSIONS 

Your questions are exactly what I like to see on this BBS. Believe me, improving DERIVE's knowl- 
edge of mathematics is much more interesting than user-interface programming. 

Your trig example is a much easier problem than the fractional power example in your previous mes- 
sage. To simplify the expression to 1, the product of trig factors needs to be collected into a single 
term. This is done by using the Manage Trig command to select "Collect” trig mode before simplify- 
ing your expression. 

Aloha, A1 Rich, Soft Warehouse, Inc. 

Message 3107: From HARALD LANG to SOFT WAREHOUSE about #3100 / SIMPLIFYING 
TRIG EXPRESSIONS 

Hmm, you caught me there. But now I am impressed. If I understand correctly, when "collect" is on, 
she (DERIVE) will rewrite COS and SIN to complex exponentials, perform the multiplication and go 
back to SIN and COS. I won't bore you with the trivial details, but she will end up simplifying the 
sum (or an equivalent sum) 

SIN(tt/18) - COS(tt/9) + COS(2tt/ 9) (1) which is =0. 

But how does she know that??? In order to outsmart her, I must give her a somewhat different prob- 
lem: simplify 

SIN(tt/7) - SIN(10tt/21) + COS(13tt/ 42) (2) which is also =0 (challenge ...) 

This one beats her, I have tried all 9 combinations of options. Note that I am not complaining on 
DERIVE, on the contrary. I will try later to sort out what I am getting at. - Harald 

Message 3109: From HARALD LANG to PUBLIC about INTERACTING WITH DERIVE 

Expression (2) in my message #3107 was SIN(7 t/ 7) - SIN(107i/21) + COS(137i/42) (2) 

This is = 0, which DERIVE can't simplify it to. I know it is equal to 0 for a peculiar reason: I looked 
at the 3 -degree equation 

x A 3 - (3/4)*x + (1/4)*SIN(37 t/ 7) = 0 (3) and it is easy to see, 

using the identities SIN(3x) = 3*SIN(x) - 4*SIN(x) A 3 and COS(3x) = 4*COS(x) A 3 - 3*COS(x), that 
this equation has the following roots: 

xl = SIN(tt/7), x2 = COS(13tt/42), x3 = -SIN(10tt/21) 

(this is also essentially the solution DERIVE gives with Manage Trig ‘Direction’^ Auto, 
’Toward’=Sines.) But since the sum of the roots of a 3-degree equation (with x A 3 -coefficient =1) 
equals the negative of the x A 2-coefficient, i.e., =0, we get my identity (2) = 0. 

Here is another way to look at it: 

consider the expression SIN(x) - SIN(x + 7i/3) + COS(x + n/6) (4) 

This is quite simple to simplify to =0, isn’t it? We just expand the second and third terms using the 
usual trig formulas. Also, DERIVE HAS NO PROBLEMS TO SIMPLIFY IT TO =0 (‘direc- 
tion’=Expand, 4 towards ’= Auto)! But if we substitute nil for x, and then simplify with DERIVE, she 
fails! Indeed, we get (2)! So we have the apparent paradox that DERIVE can simplify the more 
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general expression (4), but not the special case (2). But surely, this is true also for us humans (at least 
we find (2) more difficult than (4), don't we?) Of course, there is a similar story behind my earlier 
‘fractional power’ example. 

When I analyse a model or a problem, I usually have some information that can not be formalised: 

I know what interpretation I have in mind, where the problem come from, why I made a certain speci- 
fication, etc. There is no way (in general) that I could give all information to the programme (in a 
source code, for instance) and then ‘turn the crank’, and hope to get a ‘solution’ relevant for the ap- 
plication. Applied mathematics just isn't that way. My examples above are of course artificial, but 
you see the point: there is no way that DERIVE could know that I got the three terms in (2) as the 
roots of the equation (3). Nor could she know that the relevant way to interpret the seemingly random 
numbers 1 07 t/ 2 1 and 137i/42 is nil + 7i/3 and nil + n/6. But I know, since I put these numbers there. 
The punchline is that I need an interactive programme: I am the master, DERIVE is my clever assis- 
tant ("DERIVE — a Mathematical Assistant"; very well put!), and I must have the option in each step 
to help her and to push her the way I want. 

[As an aside: I think SWH is making a mistake in marketing DERIVE: you focus so much on it as a 
tool for education that it gives the impression that if you want something for research in applied 
mathematics, you should look for something else. But this is totally wrong: DERIVE must be (I be- 
lieve, I have no experience) a wonderful tool for education, since one can explore mathematics with it, 
not just ‘turn a crank’, but this is also precisely why it is so useful in research! (my experience is from 
economics, but isn't it generally true?) Don’t pretend that with DERIVE, you just formulate the prob- 
lem and lean back while she does all the job for you, because that is not true! Some people get disap- 
pointed and give up (I know examples of this.) It can’t be true either, this is an inherent property of 
mathematics, NOT a limitation to this very software; thus nothing to be ashamed of. And the right 
way to adopt to this reality is to make the programme interactive, as DERIVE.] 

But for me to work in an efficient way interactively with DERIVE, I must know how she "thinks", i.e., 
what she does (tries) when she simplifies, calculates an 'antiderivative', solves an equation, etc. I 
would like to hear on this BBS what experiences other users have on this, and how they manage to 
push DERIVE to perform what they want in certain situations, etc. Furthermore, SWH: I would defi- 
nitely be willing to pay some for an additional documentation on how DERIVE "thinks", I think I 
would gain at least as much efficiency from that as from an upgrade to XM. What do you say, SWH, 
any hope for that? 

Finally, here is an example I think DERIVE should manage to simplify to =0: 

(sqrt(5)+2) A (l/3) + (8*sqrt(5)+16) A (l/3) - (27* sqrt(5)+54) A (l/3) 

I have tried different option settings, but DERIVE fails. Again, one can help her by substituting x for 
sqrt(5); then she gets =0. Again it is easier for her to simplify a more general expression than a special 
case! SWH: I hope you realise how useful it would be if one could anticipate such things. -Harald 

Message 3110: From SOFT WAREHOUSE to HARALD LANG about #3107 / SIMPLIFYING 
TRIG EXPRESSIONS 

In your message #3107 you wondered how DERIVE is able to simplify 
SIN(te/18) - COS(te/9) + COS(2tt/ 9) to 0. 
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The more fundamental question is how DERIVE simplifies SIN(7i/18) + COS(27i/9) to COS(7i/9). 

Since COS(x) = SIN(x + tt/ 2) for all x, internally DERIVE stores the above expression as 
SIN(7 t/18) + SIN(137t/18) 

Then DERIVE applies the trig identity SIN(x) + SIN(y) - 2 * SIN((x + y)/2) * COS((x - y)/2) 
to the above expression. This results in the product 2 * SIN(77i/18) * COS(-7i/3) 

The above SIN factor simplifies to COS(7i/9) and the COS factor simplifies to 1/2, which then conven- 
iently cancels the 2 factor. Your done! No magic required. 

The real trick to doing computer algebra is to know which identities to apply and when to apply them. 
For example, the above trig identity should not be applied if the result is more complicated than the 
original expression. 

A bug prevented DERIVE from simplifying your second example using the same algorithm. The 
problem will be fixed in the next release. Examples like these are invaluable for finding deficiencies 
and improving DERIVE. Keep unf coming! 

Note that DERIVE does not transform SIN and COS to complex exponentials when in trig "Collect" 
mode. Rather it applies the identities listed on the bottom of page 98, Section 6.3, of the DERIVE 
User Manual. 

Aloha, A1 Rich, Soft Warehouse, Inc. 

Message 3136: From JERRY GLYNN to PUBLIC about ALGEBRAIC EXPANSION 

We think of expanding (x + 1) A 3 as algebra and when we do Taylor expansions we think we're doing 
calculus and we generate polynomials for sin(x) and exp(x) and others. Try simplifying 
Taylor((x+l) A 3,x,0,3). Looks like we can do algebra with calculus tools. It might also make better 
sense for students to start with a familiar expansion when first doing taylor. 

(x+l) A 3=aO+al*x+a2*x A 2+a3*x A 3 and differentiate the whole thing and simplify and substitute 0 for 
x and solve to work out the coefficients. 

Message 3170: From SOFT WAREHOUSE to PUBLIC about ROOTS OF UNITY 

The following definition returns a vector of the nth roots of unity ordered counter-clockwise in the 
complex plain beginning with the real root: 

l - 77 - 2 - k / ri 

NTH_R00TS_0 F_U N ITY ( n ) := VECTOR (e , k, 0, n - 1) 

For example, to compute the 5th roots of unity enter and simplify the expression 
NTH ROOTS OF UNITY (5) giving 

NTH_RQGT5_0 F_UNITY(5 ) 

JS 1 L-JQ-JS + 10 ) Js 1 L-JO-O - 2 -A) JS 1 L-JO-Q - 2 -^ 5 ) 

4 4 4 *4 4 4 '44 4 

ft 1 + 10 ) 1 


4 


4 


4 
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Message 3172: From SOFT WAREHOUSE to JERRY GLYNN about MULTIPLE EQUATION 
SOLVING 

The following definition is in response to Jerry Glynn's request for a function that individually solves 
each equation in a list of equations and returns the result as a list of lists of roots: 

SOLVE_LIST (v, x) : =VECTOR (SOLVE (ELEMENT (v, n_) , x) , n_, 1 , DIMENSION (v) ) 

Then for example, simplifying SOLVE_LIST([x-l,x A 2-l,x A 3-l,x A 4-l],x) results in the list of lists: 

[[x=l], [x=l,x=-l], 

[x=l , x=-l/ 2-SQRT ( 3 ) *#i/2 , x=-l/ 2+SQRT ( 3 ) *#i/2] , 

[x=l , x=-l , x=-#i , x=#i] ] 

SOLVE from 1994 has changed its behaviour from 1994 up to now. As you know there is a 
difference between SOLVE and SOLUTIONS. We can also write the VECTOR-expression in 
a more compact form: 


#1: 

equs := Lx - 1, x 

#2: 

SOLVE_LIST (v , x) 

#3: 

SOLVE_LIST(equs) 

#4: 

x = 1, x = -1 V 


[ 2 3 4 1 

Lx - 1, x - 1, x - 1, x - lj 


n_ 


1 73 -i. 1 73 -i 

— - v x = - — + v x = 1, x 

2 2 2 2 


-L V 


X = LVX = -lVX = l 


SOLUTIONS_LIST(v, x) VECTOR ( SOLUTIONS (v , x), n_, 1, DIMENSION(v)) 

#5: n_ 


#6: S0LUTI0NS_LIST (equs) 


#7: 



1 73 -L 1 73 -L 


[i], [i, -i], 

+ P 

, [1, -1, L, -L] 


2 2 2 2 



#8: S0LVE_l new(v, x) := VECT0R(S0LVE(v_, x), v_, v) 

#9: SOLVE_L_new(equs) 


# 10 : 


1 73 • i. 1 73 • l 

x = 1, x = -1 v x = 1, x = - — - v x = - — + v x = 1, x = -L v 

2 2 2 2 


X = LVX = -lVX = l 


#11: SOLS_l new(v , x) := VECTOR(SOLUTIONS(v_, x), v_, v) 

#12: SOLS_l new(equs, x) 

1 73 -L 1 73 • £. 

1, + , 

2 2 2 2 


HI, 11, -11, 


> [1, -1, i, -i] 


#13: 
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During the last years I received a lot of letters and I had a lot of talks and very often they ended: 

It would be fine if we could in DERIVE. 

I repeat some of your ideas: 

Functions like real(x), integer(x), ; GCD and LCM for polynomials; the commands of the menu 

should be “official" DERIVE-commands, eg. FACTOR(expression),....; the parameters of the Options 
should be implemented in commands; let's have a command to declare a variable's domain; (Dr. Kun- 
zer, Bruneck, Italy in his last letter). 

We want to create our own menus in German (or any other language) and without some commands - 
our students shouldn't be able to use soLve in the beginning; 

I wished I had an easy way to create own DEMO-files; 

Couldn't it be possible to have a RENUMBER command to avoid a line number chaos? 

When I want to load a .MTH - file I often have to use MANAGE>EXECUTE>DIR to remember the file's 
name .... 

I need another software to insert the scale 

I cannot explain my students why x = 3 in DERIVE is a horizontal line!! 
and, and, and 

Now please look at the Hardcopy: that is DERIVE 3.0: 

I think you will see many of your ideas realized. I am sure that the DERIVE USER GROUP as a whole 
took an important part in improving DERIVE. And I'm also sure that a product like DERIVE never will 
be perfect - and I hope so - because we and the ..Fathers of DERIVE" need the challenge to go on. 


#1: 

y = 2 


ttZ: 

x — —1 


#3: 

InputMode 

:= 1 


r 3 

n 

#4: 

APPROX — 

— , 


L J2 

#5: 

6.664 



2 

3 

#6: 

x - x y + y 

#7: 

m 


#8: 

u : t Rea 1 

[-1 


NUMERATOR 


#9: 


Z 



L x 



COMMAND: ^ultipl izieren Zerlegen Gehezu gleichUng Hilfe erSetzen 

Berechnen Ubertragen Fenster Loschen eMde zeiChnen elnstel lim gen 
Enter option Deriue XM 

User Free:100x Algebra 


I've changed the ALGEBRA-Menu for DERIVE novices (in German). ‘erSetzen’ is a shortcut for Man- 
age Substitute, some command the students don't need have been omitted (= hidden) (Calculus, 

Manage ). Expression #6 is an implicit function with, you can see the respective plot in the Plot 

Window. 

See also on page 17 a short DERIVE 3.0 Listing. Look there at the effect of the new Annotate Option 

If you would like to have more information about DERIVE's new features, then call or write. I'll send a 
copy of the Version 3 Release Notes. 
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Dear colleague 

I am a lecturer at the Faculty of Education. I have a Ph.D. in Mathematics and a Ph.D. in Computer 
Science. My works are mainly about applications of computers to Mathematics and about Teaching 
Mathematics with computers 

Eugenio Roanes 

I met Eugenio in Plymouth and we had lots of fun together. Eugenio is an excellent 
tennis player and a great milk drinker and he is an expert in railways. In the Dart- 
moor Nature Park he gave a lecture about the railway there when we walked along 
an old railway line. 

Some of Eugenio 's functions are still implemented in later DERIVE versions. Nev- 
ertheless they are a fine example of sophisticated programming with DERIVE es- 
pecially using its recursive capabilities. Muchas gracias Eugenio! 


DEVELOPMENT OF SOME ALGEBRAIC 
OPERATIONS ON POLYNOMIALS IN DERIVE 

Eugenio Roanes, Madrid, Spain 

Introduction. 

Derive is nowadays a common assistant for teaching mathematics both at high schools and in the first 
years of university courses. We understand that it is intended for use as an algebraic calculator (it sup- 
ports exact arithmetic) with graphical capabilities. 

We think that, because of this orientation, it is somewhat lacking, especially in polynomial handling. 
We aim to show in this paper how we have tried to overcome this, using the possibility of defining 
new functions of the 2.xx versions. 

We know that, as we cannot access internal representations of polynomials, there will be no quick way 
to obtain, for instance, a coefficient. However, there are times when these functions are really needed. 
Implementations of the functions of §2.1 and §2.2 are included in the file MISC.MTH, but they use 
differentiation. 

1 Some new integer functions. 

1.1 Division of integers 

DERIVE possesses a division function in exact arithmetic, but doesn't include a division of integers. 
A possible way of calculating the remainder of dividing two positive integers would be to use the fol- 
lowing: 


Algorithm: 

If a < b: RESTO (a,b) = a 

In other cases: RESTO (a, b) = RESTO (a-b, b) 

Implementation : 

resto(a, b) := 

If a < b 
#1: a 

res to (a - b, b) 

#2: [resto(3S4 , 35), resto(35, 384)]= [34, 35] 
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Version 2.52 includes the function FLOOR(number) 
should be better to define: 


#3: 


resto2(a, b) := a - b- FLOOR 


\ b 


- version 2.03 doesn't. Using the former, it 


#4: [re5to2(3S4 , 3B) , restD2(3B, 384)]= [34, 3B] 


1.2 Greatest Common Divisor. 

The above makes it possible to define a function GCD in N using Euclid's algorithm. We have done it 
this way: 

Algorithm: 

MCD(a,b) decides which of the entries is greater and sends to 
AUX (greater { a, b }, smaller { a, b } ) . 

AUX(a,b) is Euclid's Algorithm: the division of integers is re- 
peated, and at every step the divisor becomes the dividend and 
the remainder becomes the divisor. When the remainder is 0, the 
previous remainder is given as output. 

Implementation : 

MCD_AUX(c , d) := 

If d = 0 
#B : c 

MCD_AUX(d , restoCc, d)) 

MCD(a, b) := 

If a > b 

#6: MCD_AUX(a, b) 

MCD_AUX(b, a) 

#7: MCD(324B6, 2SBS8) = 4 

We do have now function GCD(a.b) - which is working for numbers only! ! 

#8: GCD(324B6, 2SBSS) = 4 


[t 1 gebralca rdotherlprgn I ole 1 ean 

upT 1 

" {res to (a - b , b), else * resto( -* ’ 

Done 

■ resto(3S4 , 35) 

34 

■ a - b ■ f loorj^J + resto2{a , b) 

Done 

■ resto2(384 ? 28) 

2Q 

■ £ mod(334 , 35) mod(334 , 28 )} £ 34 

20} 

<mod<384.35}.mod<384.28>M 

DNLiE RAD liUTD FUNC E/30 



And there is a function on the handheld, too 


2 Some new functions for univariate polynomials. 

Unfortunately, DERIVE only allows the substitution of a variable by a certain value in a function of a 
polynomial if the user is in the interactive mode. So, we can only find the possibility of defining a 
substituting function: SUB (polynomial, variable, value) as the limit of the function “polynomial 46 
when ‘Variable 44 tends towards “value 44 : 

sub(p, v, a) := lim p 
#9: v-»a 

In the meanwhile function SUBST is available and we don't need sub(p,v,a) any longer. 
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2.1 Degree of a polynomial. 

The degree 0 term of a polynomial p(x) is p(0). Then p(x) - p(0) is divisible by x. The times we can 
repeat this process until we obtain 0 is the degree of the original polynomial. 

Implementation : 

p - SUBST(p, V, 0) 

#10: POLJUJXCp, v) := 


GRAD0(p, v) := 

If p = 0 
# 11 : -1 

1 + GRAD0(P0L_AUX(p, v) , v) 
1 + GRAD0(P0L_AUX(p, v) , v) 


3 2 4 

#12: GRADO (3-x + 2-x - 1)- 3-x - 3-x + 2-x 


2.2 Coefficient of the term of degree n. 

Algorithm: 

Po (x) = p (x) 

Pn(x) = (p n -i(x) - p n -i(0))/x 
COEFF (p (x) , n) = p n (0) 

Implementation : 

3 2 4 

#13: (3-x + 2-x - 1)- 3-x - 3-x + 2-x 


S 5 

3-x 7 6 9-x 4 

#14 : - + 6 ■ x - x - + 7 ■ x 




P0L_AUX2(p, v, n) := 

If n = 0 
#15 : p 

(P0L_AUX2(p, v, n - 1) - SUBST(P0L_AUX2Cp, v, n - 1) , v 
#16: COEFFCp, v, n) := SUBST(P0L_AUX2(p, v, n) , v, 0) 


#17: 

COEFF 

3 

(3-x + 2-> 

E - !)■ 

r b 'i 

2 4 x 

3-x - 3-x + 2-x - 

, x, 5 





L 2 J 

J 


LCQEFFCp, v) := COEFFCp, 


#15: 


LCOEFF 
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p - REST0_P0L( P| q, v) 

#27: C0C_P0L(p, q, v) := 

q 

To calculate the pseudoremainder we shall multiply the dividend by the leader coefficient of the divi- 
sor to the power of the difference of degrees plus one. This way we can ensure that if the polynomials 
belong to Z(x), no rational coefficient will appear in the calculations. 

Implementation : 

1 + GRADOCp, v) - GRADOCq, v) 

#28: SREST0_P0L( P| q , v) := REST0_P0L (LCQEFF (q , v) -p, q, v) 

7 5 4 2 3 2 

#29: CQC_PQL(3-x - x + 3-x - x - 2, 4-x -3-x + x + 1, x) 

4 3 2 

768 -x + 576 -x - 16-x + 420-x + 175 

#30: 

1024 


7 5 4 2 3 2 

#31: RE5TO_P0L(3-x - x + 3-x - x - 2, 4-x - 3-x + x + 1, x) 

2 

903-x + 595 -x + 2223 

#32: 

1024 


7 5 4 2 3 2 

#33: SREST0_P0L(3-x - x + 3-x - x - 2, 4-x - 3-x + x + 1, x) 

2 

#34: - 903-x - 595-x - 2223 

(The computing times were about 10 seconds in 1994. Compare with computing times of today!!) 


2.5 GCD of two polynomials 

Euclid's Algorithm also works for calculating the GCD of two polynomials. The only changes which 
have to be made are: 

• Order the polynomials in the first step according to their degree, 

• Substitute the function: remainder of the division of integers by the function: remainder of the divi- 
sion of polynomials. 

One disadvantage of this algorithm is that the coefficients may increase greatly in size. 

Implementation : 

MCD_PQL_AUX(c , d, v) 

If d = 0 
#35: c 

MCD_P0L_AUX(d , RESTQ_PQL(c , d, v) , v) 

MCD_P0L_AUX(d , RESTQ_PQL(c , d, v) , v) 

MCD_P0L(a, b, v) 

If GRADQCa, v) > GRADQ(b, v) 

#36: MCD_P0L_AUX(a, b, v) 

MCD_P0L_AUX(a, b, v) 

32 32 144500 

#37: MCD_P0L(x - S-x + 21-x - 18 , x + 9-x + 27-x + 27, x) = 

3249 

5432 432 2 

#38: MCD_P0L(x -5-x +4-x +x - 5-x t 4, x -5-x + 5-x -5-x+4,x) = 2-x -10-x+S 

2 

43 2 3 2 34- (x - 5-x + 6) 

MCD_P0L(x - 7-x + 18-x - 22-x + 12, 3-x - 13-x + S-x + 12, x) = 


#39: 


9 
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The pseudoremainder could be used for polynomials of Z(x) (instead of the remainder). In this way all 
the polynomials in every step of the algorithm should have integer coefficients: 

Implementation : 

MCD_POL_AUX_Z(c , d, v) : = 

If d = 0 
#40: c 

MCD_P0L_AUX_Z(d , SRESTO_P0L(c , d, v) , v) 

MCD_P0L_AUX_Z(d , SRESTO_P0L(c , d, v) , v) 

MCD_P0L_Z(a, b, v) := 

If GRADOCa, v) > GRAD0(b, v) 

#41: MCD_P0L_AUX_Z (a , b, v) 

MCD_P0L_AUX_Z (a , b, v) 

4 3 2 3 2 2 

#42: MCD_P0L_Z(x - 7x + IB-x - 22-x + 12, 3-x - 13-x + S-x + 12, x) = 34- (x - B-x + 6) 


Bibliography 
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Now with DERIVE 3.0 and higher ... 





3 

r b ^ 

2 4 x 


#43: 

POL 

YJIOEFF 

(3-x + 2-x - 1). 

V. 

3-x - 3-x + 2-x 

L 2 J 

i x, B 

J 


4 3 2 3 2 2 

#44: MCD_PQL_Z (x - 7-x + IB-x - 22-x + 12, 3-x - 13-x + S-x + 12, x) = 34- (x - B-x + 6) 


4 3 2 3 2 2 

#43: PQLY_GCD(x - 7-x + IS-x - 22-x + 12, 3-x - 13-x + S-x + 12) = x - B-x + 6 

4 3 2 

7 B 4 2 3 2 3 x 9x x 10B-x 17B 

#46: qU0TIENT(3-x - x + 3-x - x - 2, 4-x - 3-x + x + 1) = + + + 

4 16 64 2B6 1024 


7 B 4 2 32 903 -x B9B-x 2223 

#47: REMAINDER(3-x - x + 3-x - x - 2, 4-x - 3-x + x + 1) = 

1024 1024 1024 

Now let me try to implement Eugenio's procedures on the Tl-handheld devices. As you can see in 
the following, they works without any problems. It is easy to repeat the recursive functions or to re- 
write them using while-loops. First I worked with the TI-92 / Voyage 200. After installing the new Tl- 
Nspire Computer Link Software I defined the functions and repeated the whole procedure also on 
the Nspire. It is a pity that it is not able to transmit the TI-92/V 200 functions to the Nspire although 
the CAS-machine is more or less quite the same. See first the results (the first screens from the V 
200 ): 
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Here are the screen shots from the Nspire. You can download the TI-92/V200 and the Nspire files as 
well. The PC-Nspire works identically. 
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RAD AUTO REAL 




RAD AUTO REAL 


irno 


RAD AUTO REAL 





Compare ruf_aux and ruf_aux2: the while-loop versus the recursive function. Which way to program 
do you prefer? 










Both programming styles lead to the correct result, restopol on the Nspire and on the Voyage 200: 


RAD AUTO REAL 


Define restop 


or,vJ=Func 

Local gddo,gdor,?ddo 
gra d o {ddo; v ) -* gddo 
gra d o {dor, v ) -* gdor 
lcoeff(^o,v)-^^^o ► 
Ic o eff {dor, v ) -* Idor 


rFlTmnY YF::tYfHtY FF YFfirY 'I 

| t j— |Contr o 1 1 1 /Q |0ar I f i nd... |Mode | | 

: restopoKddo, dor, oi 

: when tgrado*; ddo, oXgradotdor, , ddo, rest 
opo K ddo-dor* 1 coef f t ddo , u X U coef f £ dor , 
u ? } £ grade £ ddo, o>-grado£dor, > , dor, o 


when sddo < sdo r, ddo 


DHL1E 


FiliD hUTD 


FUNC 
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The following comments are obsolete, because of changing the output of the SOLVE- 
command and introducing SOLUTIONS. But it might be of historical interest. Josef 


Comments on the RHS-function 

J.Bohm, Wiirmla, Wiirmla, Austria 

V.Neurath asks in his letter from 15.06.94 to explain the “undocumented" RHS- and LHS-function I 
used in a file which I've sent him. RHS(relation) and LHS(relation) returns the righthand side or the 
lefthand side of a relation. The two functions are available since version 2.57. In the beginning I 
didn't know where to use these functions but another question of V.Neurath and a letter from Dr. 
Kuenzer brought some ideas. I wanted to write a “DERIVE-program" to investigate a curve. 

Many thanks for Dr.Kuenzer's ideas (Symmetry, ....) 

See first an example to illustrate the RHS-function: (This is DERIVE for DOS from 1994!) 

1: "Solution of equation equ for variable u: M 

2: SOL (equ, u) := SOLVE (equ, u) 

2 

3: SOL (x - 4 x + 5, x) 

4: [x = 2-i, x = 2 + i] 

5: "Using RHS we select the value of the 2nd solution:" 

6: RHS (ELEMENT ( [x = 2 - 1, x » 2 + i] , 2 ) ) 

7: 2 + i 

8: "See the vector of the solutions (without x =) ! Now I can use the solutions" 

9: VEC_SOL (equ, u) := VECTOR (RHS (ELEMENT ( SOLVE (equ, u) , i)), i, 1, 

DIMENSION (SOLVE (equ, u) ) ) 

2 

10: VEC_SOL (x - 4 x + 5, x) 

11: [2 - i, 2 + i] 

12: "I want to substitute the solutions of equ for variable v in expression expr: " 

SUBST (equ, u, expr, v) := VECTOR ( lim expr, i, 1, 

13: v->RHS (ELEMENT (SOLVE (equ, u), i)) 

DIMENSION (SOLVE (equ, u) ) ) 

2 2 

14: SUBST (x - 4 x + 5, x, 2 a - a + 2, a) 

15: [6 - 7 i, 6 + 7 i] 

The appearance in modern DERIVE: 

2 

#1: SOLVE (x - 4-x + S, x) = (x = 2 - i. v x = 2 + 0 

2 

#2: RHS(SQLVE (x - 4-x + S , x)) = (x = 2 +■ l) 

2 

#3: RHS(RHS(SOLVE (x - 4-x + S, x))) =2 + l 

2 

#4: SOLUTIONSCx - 4-x t 5 „ x) = [2 + l , 2 - l] 

2 

#B : VECTOR (2- a - a + 2, a, [2 + l, 2 - LJ) = [6 + 7 - l , 6 - 7-1] 

#6: substi (equ , u, expr, v) := VECT0R(expr, v, S0LUTI0NS(equ , u)) 

2 2 

substi'Cx - 4-x + B, x, 2-a - a + 2, a) = [6 + 7 - 1 , 6 - 7-i.] 


# 7 : 
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One of the many reasons I like DERIVE - it is not the most important one - is that it is a pow- 
erful tool for my work as schoolteacher. V.Neurath's and Dr. Kuenzer's letter made me curi- 
ous if it could be possible to write a sort of program for “discussing a curve”. 

Some hours later I saw the product of my efforts and I was satisfied. Here you can see the 
listing. I think, it doesn't need any further explication. The function has to be edited within 
brackets, when it is a rational function then enter f := [numerator, denominator], because I 
don't use version 3.xx. 

This was the leading paragraph in DNL#15. I reprint the first lines of codes. Then I will proceed with a 
recent version of the “program” which is not really a program, but a collection of functions. Really pro- 
gramming became true with version 5. Josef 

"Discussion of a curve or: how to write a ’program’ with DERIVE!" 

1: f: = [] 

2 : z : =ELEMENT (f , 1) 

3 : n:=IF (DIMENSION (f ) =1, 1 , ELEMENT (f , 2) ) 

4: [nl : =DIF (n, x, 1 ) , n2 : =DIF (n, x, 2 ) , zl : =DIF ( z , x, 1 ) , z2 : =DIF ( z , x, 2 ) , f 1 : =z/n] 

5: "Zeros" 

6: nO : =SOLVE (z,x) 

7 : zer : =APPEND ( [ ["Zeros [fl=0, ""] ] , VECTOR ( [IF (IM(RHS (ELEMENT (nO, i) ) ) , 

RHS (ELEMENT (nO, i) ) , "komplex") , IF ( IM (RHS (ELEMENT (nO, i) ) ) , 0 ,"!")], i , 1 , 
DIMENSION (nO) ) ) 


Discussion of a cun/e or 

How to write a "program" with DERIVE! 

#1: [Preen si onDi gits := 6, NotatnonDign ts := 6] 

#2: f(x) := 

#3 : [z(x) := NUMERATOR (FACTOR (f (x) , Trivial)), n(x) := DENOMINATOR (FACTOR (f (x) , Trivial))] 

f ffd Tl T ffd T2 T] 


#4: 


fl(x) := FACTOR — f (x) , Tri vi al , f2(x) := FACTOR — 


f (x) , Trivial 


LdxJ 


LdxJ 


#5: 

d 

fll :z — f(x) , 

f22 := fi? f(x) 


dx 

LdxJ 


Finding the zeros: 

#6: nO := SOLUTIONS(z(x) , x, Real) 


#7: 


APPEND 


Zeros : 
f(x) = 0 


, VECTOR ([v , 0], v, nO) 


#S: f(x) 


4 3 2 

x - 11- x + 29-x + 35-x - 150 


30 


#9 : zer 


Zeros : 


3 2 

x- (x - 11- x + 29 -x + 35) = 150 


-2 

3 

5 


0 

0 

0 
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Finding the local extremal values: 

#10: eO := SOLUTIONS (NUMERATOR (fl(x)) = 0, x, Real) 


^281 13 13 7281 


#11: eO = 5, 


#12: extr :: APPEND 


Extremals: 


f'(x) r fl(x) 


, VECTOR C[v, f(v), IF(SIGN(SUBST(f22, x, v)) = 1, MIN, MAX, 


TerrPt)], v, eO) 


Extremals: 


f'Cx) = 


3 2 

4-x - 33 ■ x + 58- x + 35 


#13: extr = 


7281 13 


843-7281 38939 


13 7281 


5120 15360 


843-7281 38939 


Finding the points of inflection: 


5120 15360 


#14: wO := SOLUTIONS (NUMERATOR (f 2 (x)) = 0, x, Real) 


#15 : wO = 


7393 11 11 7393 


#16: inf 1 := APPEND 


Infl . points: 

f"Cx) = f2(x) 


VECTOR (Tv, f(v), lim fll 


j, v, wO) 


Infl . points: 


f"Cx) = 


6-x - 33 -x + 29 


#17: infl - 7393 11 5-7393 12457 15 131-7393 


8640 16 


11 7393 5-7393 12457 131-7393 15 


L 4 12 64 8640 2160 16 J 

Finding the poles (vertical asymptotes): 

#13: pO := SOLUTIONSCnCx), x) 

#19: pol := APPEND([[Poles:]], IF(DIMENSIGN(pO) = 0, [[none]], [pO] 1 )) 

Finding the endbehaviour: 

#20: as_ := QUOTIENTtz(x) , n(x)) 

#21: as := IF(as_ - f(x) = 0, none, as_, as_) 

#22: ass := APPENDS [Asymptote: ]] , [[as]]) 

Determining the symmetry 


#23: sym := APPEND([ [Symmetry : ]] , IF(f(x) - f (-x) , [[axial symm.]], IF(f(x) + 


f(-x), [[centr.symm. ]] , [[no symm.]], [[no symm.]]), IF(f(x) +f(-x), 
[[centr . symm. ]] , [[no symm.]], [[no symm.]]))) 
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Finding the limits forx -»±«: 

"linns APPEND([[Li'nn' ts : , ]], 


— > co. "h'm f(x 


*]]■ [[■ 


extr 

infl 

#25: discuss := pol 


#26: discuss 


x -O - 11 -x + 29-x + 35) = 150 


Extremals : 


4 ■ x - 33 ■ x + 5 8 ■ x + 35 


1 A - 

30 


5 

0 


^281 13 

843-^281 

38939 

8 8 

o 

r-j 

■ — i 
LTi 

15360 

13 ^281 

843-^281 

38939 


^120 15360 


Inf 1 . points : 


f rr 00 = 


6 ■ x - 33 ■ x + 29 
15 


^393 

11 

5 -^393 

12457 

15 

12 

4 

64 

8640 

16 

11 

^393 

5 -^393 

r- 

LTi 

r-j 

■ — i 

1 

131 ^3! 

4 

12 

64 

8640 

2160 


L none J 
Limits : 
x — > 

x — > & 

Asymptote: 
none 

Symmetry: 

. no symm . 


You can approximate #27 and then you will receive decimal numbers. Use the results to plot 
the graph including all critical points. 
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Graph of the function from above followed by 
a 2 nd example: 


5 x 

fO) := 

#33: x 2 


x - 1 


#34: discuss 


#56: 



Zeros : 
x 5 


2 

x - 1 


1,11805 

-1,11805 


Extremals : 


4 2 

4 ■ x - 11 -x +5 


f r 00 = 

2 2 2 
x -O - 1 ) 


1.47492 

2.15517 

MAX 

-1.47492 

-2.15517 

MIN 

0.758029 

8.578 

MIN 

-0.758029 

-8.578 

MAX 


Inf 1 . points : 

f rr W = 

1.88165 

-1.88165 


slope 


6 4 2 

2 ■ (4 ■ x - 18-x + 15 -x - 5 ) 

5 2 5 

x -O - 1) 


1.91662 
-1.91662 
Poles : 

0 
1 

-1 
Limits : 
x — > ,*1 0 

x — > -,*1 0 

Asymptote : 

0 

Symmetry : 
centr . symm . 


-0.708727 

-0.708727 


Try 


- x/2 2 

f(x) := e ■ (x - 4) 


or any polynomial or rational function which are usually treated in school. 


p 26 


Josef Bohm: 


Discussion of a curve 


D-N-L#15 


Graph of the function from above and a final 
example: 

3 

H + H - 1 

#45: fO) := 

1 - 5 -x 

#46: discuss 
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Aufgaben aus der Elektrotechnik - nicht nur 
fur den Unterricht in beruflichen Schulformen 2 


H. Scheuermann u. N.Wild, Hofheim, Germany 


5. Netzwerkberechnung - einmal anders 
(Calculation of networks - another way) 8 

Aufgabenstellung: 

Berechnung eines beliebigen elektrischen Wider- 
standsnetzwerkes mit Hilfe eines Naherungs- 
verfahrens. Ein triviales Beispiel soil den Algorithmus 
erlautern und das iterative Vorgehen demonstrieren: 

Problem: Calculation of an arbitrary network of resis- 
tors using an approximation method. A trivial example 
shall explain the algorithm and demonstrate the iterative 
process: 



Fur die Schaltung in Abb. 8 gelten die folgenden Werte: 


U = 10V; 
R 2 = 20Q; 


R-, = 10Q 
R 3 = 300 


Maschengleichungen (Mesh-equations): 

(I) hRi + l 2 R 2 - U = 0 

(II) l 2 R 2 + l 3 R 3 = 0 


Abb 8: Lineares Widerstandswerk I 

Linear network of resistors 


lm Gegensatz zu den Losungsverfahren wie z.B. Maschenstrom verfahren, Superposition nach Helm- 
holtz, Ersatzspannungs- und Stromquelle usw., wird bei diesem Verfahren zur Netzberechnung nur 
ein (beliebiger) Pfad vom positiven bis zum negativen Anschluftpol der Spannungsquelle ausgewahlt. 
Alle Strome dieses Weges erhalten einen beliebigen Startwert ungleich 0, zB 1, alle anderen Strome 
den Wert 0. 

ii = i 2 = i; i 3 = 0 

Mit diesen Werten sind die Maschengleichungen (I) und (II) naturlich noch nicht erfullt. Deshalb wer- 
den die Gleichungen durch einen additiven Korrekturwert Ai erganzt. 

We choose one - optional - path from the positive to the negative connection terminal of the supply point. All 
the currents along this way get an arbitrary starting value ^ 0, e.g. 1 , the other currents have the value 0. These 
values don't meet the equations (I) and (II). So we add a correcting value Aj in (I) and (II).) 




U - I l R l- I 2 R 2 

V R 2 


1 : = I 1+ A 1 

1 neu 1 1 

T : = I. + A 

2 neu 2 1 


(without DERIVE ! ! ! ! , editor) 


Die zweite Masche wird nun entsprechend korrigiert: 


I. + A R 
2 2/2 


r 3- A 2j R 3 


l 3 R 3 l 2 R 2 
R 2+ r 3 


E : = E 
2 neu 2neu 


I. - A. 

3 2 


Die Idee zu diesem Verfahren wurde von Klingen [3] veroffentlich. Klingen [3] published the idea of this method. 
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It I, It, 

' 10 20 o 

#1: r_mat : = 

. 0 20 -30 

"11 0 0 

#2 : e_m at : = 

. 01-10 

#3: i_vec := [1, 1, 0, 1] 

Die Korrekturwerte A fur die i-te Zeile liefert nach Gleichung (*) der Ausdruck: 

R (i) • I 

A. = - — - — (the correcting values for line i) 

R (i) E (0 

Damit la(Jt sich mit DERIVE eine Funktion der Korrekturwerte angeben: 

r_mat *i_ 
n 

#4: S(n_, n) := 

r_mat -e_mat 
n n 

Die gesamte Berechnung einer Naherung kann durch die folgende Beziehung beschrieben werden, 
wobei i die Werte von 1 und 2 (Anzahl der Maschen) durchlauft: 

Inew-IoM " A . E (j) (1=1,2) 

Die von DERIVE bereitgestellte Schleifenkonstruktion ITERATE bietet sich zur Berechnung des Aus- 
drucks an, kann aber nicht verwendet werden, da l new durch zwei Parameter l oid und i bestimmt wird, 
ITERATE aber nur einen Parameter ubergeben kann. Da sich jede Iteration auch durch eine Rekursi- 
on beschreiben laftt, kann auf die folgende Programmiertechnik zuruckgegriffen werden: 



In DERIVE lautet diese Konstruktion: 

I O' _ i n) := 

If n < 2 

#B : I(W_ - &(1_i n) ■ e_mat|.n , n + 1) 

i_ 

Mit dem Aufruf l(i_vec,1) liefert die Funktion I dann eine Naherung. Weitere Naherungen lassen sich 
durch die ITERATES - Funktion generieren: 

(Now we have only one parameter i_, so we can use the ITERATES-function:) 

#6: ITERATES(I(i_, 1), n_, i.vec, 6) 

Hierbei gibt der erste Parameter die zu wiederholende Funktion (also l(i_, 1 )) an, der zweite das Funk- 
tionsargument (i_), der dritte den Startwert und der vierte die Anzahl der Wiederholungen. 

Die Ausfuhrung der Zeile #6 liefert (Approximating #6 gives): 


U 

(the coefficients of (I)) 
(the coefficients of (II) 
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1 

1 

0 

1 

0.333333 

0.2 

0.133333 

1 

0.422222 

0.253333 

0.16&&&& 

1 

0.445925 

0.267555 

0,17837 

1 

0.452246 

0.271348 

0.180898 

1 

0.453932 

0.272359 

0.181573 

1 

0.454381 

0.272629 

0.181752 

1 


Das Losungsverfahren ist sicherlich nicht so potent, daft es die genannten traditionellen Berech- 
nungsverfahren verdrangen konnte. Dennoch sprechen auch Grunde dafur, dieses Verfahren - insbe- 
sondere neben den klassischen - in den Unterricht zu integrieren: 

1. Ein uberraschender und fur die Schule unublicher, in der Praxis aber durchaus gebrauchlicher 
Losungsansatz (Motto: erst vereinfachen wir die Losung und machen dabei Fehler, die wir spater 
sukzessive wieder korrigieren). 

2. Wenn die Schuler schon Kenntnisse und Ubung in der Anwendung von DERIVE haben, dann stellt 
die Aufgabe ein weiteres Beispiel dar, daft Operationen der Vektoralgebra auch zur Losung von 
nicht-geometrischen Problemen verwendet werden konnen. 

This method is not so powerful to substitute traditional calculation methods. But from my point of view there are 
two reasons to integrate this method into classes: 

1 . A surprising - unusual in school but very usual in practice - approach (at first we simplify the solution risk- 
ing errors, which we then correct successively). 

2. If the students have some skill in working with DERIVE this problem can be a fine example to show that 
vector algebra operations can be used successfully to solve non geometric problems. 

With Derive 5 and higher we are able to collect the procedure into one program and overcome 
Helmut's complaints from above. The code line after the “Loop” (i :=[i>kL, ...]) is necessary. Try the 
program without this line to experience the difference. Josef 


"lin_net(r, e, i, n, its, 51, 52, k) 
Pnog 

its := [f] 
k := 0 
Loop 

-i := [ul, - 42 , 43, 1] 
its := APPENDCWts, [i]) 

If k = n 

RETURN itsitl n + 1] 

51 := - rj,l- 1 /(rjJ.- ej,l) 

41 : = 41 + Si 

42 := 42 + 51 

52 := - rj,2- 1 ej,2) 

42 := 42 + 52 

43 := 43 - 52 
k :+ 1 


~h'n_net(r_mat , 
1 

0.3333333333 
0.4222222222 
0.4459259259 
0.4522469135 
0.4539325102 
. 0.4543820027 


e_mat, i_vec , 
1 

0.2 

0.2533333333 
0.2675555555 
0. 2713481481 
0.2723595061 
0.2726292016 


6) 

0 1 
0.1333333333 1 

0 . 16S8SSSS8S 1 
0.1783703703 1 

0.1808987654 1 

0.1815730041 1 
0.1817528010 1 


And what about the handheld? 


PiSklglEra 

( F3T Y FHt y F! Y FfiT Y 1 

|Calc|0ther|PrgnI0|dean Up| | 


■ i . 000000 

1 . 000000 

G. GGOGGc 

5 1.000000' 


. 333333 

. 2GQQ00 

. 133333 

1 . 000000 


. 422222 

. 253333 

. 168889 

1. 0GGG00 


. 445926 

. 267556 

. 17837G 

1 . 000000 


. 452247 

.271348 

. 180899 

1 . 000000 


. 453933 

. 272360 

. 181573 

1. 0GGG00 


■ 4543S2 

■272629 

■ 181753 

1.000000 

lin_net <r 

_mat ,e_ 

.mat ,i_ 

uec , b 

1 DNL1S RAD AUTO 

FUNC 1 E/30 1 


(» Ifl 1 gebrafca fclotherfprgr'i I dIc 1 ean Up! 1 

r 10 20 

0 - 1 G 1 




■ 

r_Fiat 



L 0 20 

-30 G J 

1 

[10 

2 G 

G - 1 G 1 



[o 

2 G 

"30 G J 

ri i g 

01 


ri 

1 0 01 

■ 

+ e_nat 




L 0 i - 

1 0 j 


Lo 

1 "I 0 j 

■[1 1 G 

1] + i_yec 


[1 

1 G 1 ] 

1 

DNL1S 

RAD AUTO 

FUNC 1E/30 
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6. Berechnung gemischter passiver Schaltungsnetzwerke 
(Calculation of mixed passive circuit networks) 

Die Berechnung von gemischten Schaltungsnetzwerken hat i.a. zwei Ziele: 

1 . Umfangreichere Schaltungen auf einfachere Ersatzschaltungen zuruckzufuhren und 

2. bestimmte (Schein-)Widerstandwerte durch die Kombination vorhandener Bauteile zu erzielen. 

(The goals are to reduce complex circuits to simpler alternative circuits and to achieve certain impedances com- 
binating the existing components.) 

DERIVE bietet schon von der Programmstruktur her geeignete Programmelemente, die eine elegante 
Losung der oben genannten Aufgabenstellungen ermoglichen (Funktionsdeklaration, fest implemen- 
tierte Funktionen, komplexe Zahlen, einfache Programmierungen usw.) 

Im Folgenden sollen zur Vorstellung und Erlauterung des Losungsweges zuerst Gleichstrom- 
Schaltungen und im weiteren auch Wechselstromschaltungen berechnet werden. 

6.1. Berechnung von linearen Schaltungen 
(Calculation of linear circuits) 

Aufgabenstellung: 

Der Gesamtwiderstand R tot der in Abb. 9 dargestellten Schaltung soil bestimmt werden. 

Calculate the total resistance of the circuit in fig. 9 

Zur Losung der Aufgabenstellung werden die beiden aus den Kirchhoffschen Gesetzen hergeleiteten 
Gleichungen „geeignet kombinierf : 



Die Berechnung kann mit Hilfe eines programmierbaren Systems, wie z.B. Taschenrechner, PC mit 
Programmiersprachen Oder mit entsprechender mathematischer Software, durchgefuhrt werden. 
Sinnvoll ist dabei, daR> nicht zuletzt aus schreibtechnischen Grunden die Berechnungsvorschrift „in 
einer Zeile“ darstellbar ist. Naturlich kann das Berechnungsverfahren auch wie oben beschrieben 
unter Verwendung von Variablen bzw. Makro(er)setzungen formuliert werden, jedoch ist dies bezug- 
lich der Rechenzeit und des Speicherbedarfs wenig vorteilhaft. In der herkommlichen Form, also mit 
Bruchen geschrieben, wird der Ausdruck fur R tot „unhandlich“: 

l 

R = 

1 

1 

1 1 
r 3 R 4 + R 5 


1 

R, 


Ro + 








P 32 


H. Scheuermann & N. Wild: Elektrotechnik 2 


D-N-L#15 


Eine elegantere Darstellung wird durch die Definition von Rechenoperationen fur die Reihen- und 
Parallelschaltung erzielt: 

'+ ' := Reihenschaltung (series conn.) 7 T := Parallelschaltung (parallel conn.) 

Damit folgt fur das Beispiel in Abb. 9: 

Rtot = Ri // ( R2 + ( R3 II ( R4 + Rs)) Infix - Notation 

Diese Schreibweise entspricht der sogenannten Infix-Notation einer zweistelligen Rechenoperation 
der Algebra, d.h. der Operator steht zwischen den Operanden. Daneben gibt es noch die Prefix- und 
die Postfix-Notationen, auch polnische (PN) bzw. umgekehrte polnische Notation (UPN) genannt. 9 

Die Schaltung in Abb. 9 kann dann wie folgt in Prefix-Darstellung geschrieben werden: 

Rtot = //; + ;//; + R5 ; R4 ; R3 ; R2 ; Ri Prefix - Notation 

In dem benutzten programmierbaren System konnen nun zwei Funktionen SER(r1,r2) := ri + r2 und 
PAR(r1 ,r2) := 1/(1/r1 + 1/r2) definiert werden. Der Gesamtwiderstand R tot ergibt sich dann aus: 

Rtot = PAR(SER(PAR(SER(R 5 , r 4 ) , R 3 ) , R 2 ) , Ri ) 

Fur die Berechnung einer Schaltung mit dem Programm DERIVE genugt die Definition der Funktionen 
SER und PAR, sowie der Term der Schaltung in Prefix - Notation. 

#1: Definition of the operations: 

#2: SER(rl, r2) := rl + r2 

1 

PAR ( rl , r2) := 

#3: 11 

+ 

rl r2 

#4: Description of the circuit: 

#B: rtot = PAR (SER (PAR (SER ( rB , r4) , r3) , r2) , rl) 

#6: general solution: 

rl-(r2-(r3 + r4 + rB) + r3- (r4 + rB)) 

#7: rtot = 

rl- ( r3 + r4 + rB) + r2-(r3 + r4 + rB) + r3-(r4 + rB) 

#S : values for rl = li"j, r2 = 2Q 

1- (2- (3 + 4 + B) + 3- (4 + B)) 

#9: rtot = 

1- (3 + 4 + B) + 2- (3 + 4 + B) + 3- (4 + B) 

17 

#10: rtot = 

21 


9 Ein einfaches Beispiel verdeutlicht die unterschiedlichen Schreibweisen: 

Infix: (3+4) * (5+6) 

Postfix: 5; 6; + ; 3; 4; +; * 

Prefix: *; +; 3; 4; +; 5; 6 

In der Prefixdarstellung wird zuerst der Ausdruck A1: +; 5; 6 und dann A2: +; 3; 4 ausgewertet. Die letzte Operation *; A2; A1 
liefert dann das Ergebnis. Diese Arbeitsweise ahnelt stark dem Aufruf von Funktionen - in diesem Fall Funktion + und Funktion * 
mit je zwei Parametern. 
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Abb 10: Bruckenschaltung (mit Dreieck-Stern-Umwandlung 


Das zweite Ziel der Berechnung von gemischten Schaltungen, einen bestimmten Widerstandswert 
aus der Kombination vorhandener Bauelemente zu erzielen, stellt im gewissen Sinn die Umkehrauf- 
gabe dar. Der Gesamtwiderstand ist gegeben und ein einzelner Widerstand ist gesucht. Mit DERIVE 
ist es kein Problem, die Gleichung fur R tot nach einem anderen Widerstand umzustellen. Der Befehl 
soLve nach R1 liefert, angewandt auf den Ausdruck #5 Oder #7, die Zeile #1 1 : 

The 2nd aim of calculating mixed circuits - to get a certain resistance value from a combination of existing com- 
ponents - is a sort of reverse task: you know the total resistance and you want to find a particular one. So soLve 
expression #5 or #7 for R1 to get expression #11: 

rtot-(r2-(r3 + r4 + rB) + r3-(r4 + rB)) 

#11: rl = 

r2-(r3 + r4 + rB) + r3-(r4 + rB - rtot) - rtot- (r4 + rB) 

Substituting for r2, r3, etc and rtot = 17/21 we obtain the original value for rl (as we have expected). 

17 

(2- (3 + 4 + B) + 3- (4 + B)) 

21 

#12 : rl = 

f 17 T 17 

2- (3 + 4 + B) + 3- 4 + B (4 + B) 

L 21 J 21 

#13 : rl = 1 

Hat das gegebene Netzwerk Elemente einer Bruckenschaltung, so kann das beschriebene Losungs- 
verfahren auch, aber nicht ohne Modifikation, angewandt werden. Durch Anwendung der Dreieck - 
Stern - Umformung kann der Gesamtwiderstand R tot der Schaltung laut Abb. 10 nach beschriebenen 
Losungsverfahren bestimmt werden. 

Mit Hilfe der definierten Funktionen SER und PAR (in den Hintergrund - als Utility file - laden) kann die 
Gleichung fur den Gesamtwiderstand der Bruckenschaltung wie folgt geschrieben werden: 
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Die Funktion TS (expr. #1) transformiert denjenigen Widerstand, der als erster Parameter in TS aufge- 
fuhrt wird. Z.B. liefert der Aufruf von TS(r2,r1 ,r3) den Widerstand R2\ 

If there are elements of a bridge circuit within the network we can use the algorithm described, too but not with- 
out a modification. Applying the Triangle - Star - modification R tot of the circuit in fig. 10 can be determined. 
Load the functions SER and PAR as Utility files and write the equation for R tot . TS transforms the first resistor 
of the parameter list, eg. TS(r2,rl,r3) leads to R2'. 

#B : SER(rl, r2) := rl + r2 

1 

PAR ( rl , r2) := 

# 6 : 11 

+ 

rl r2 

r2- r3 

#7: TS(rl, r2 , r3) := 

rl + r2 + r3 

#S : rtot = SER (PAR (SER (rl, TS(r4, r2 , rB)), SER(r3, TS(r2, r4 , rB))), TS(rB, r2 , r4)) 

rl- (r2- (r3 + r4) + r3- (r4 + rB) + r4- rB) + r2-(r3-(r4 + rB) + r4-rB) 

#9: rtot = 

rl-(r2 + r4 + rB) + r2- (r3 + rB) + r3-(r4 + rB) + r4-rB 


6.2. Berechnung von Wechselstromkreisen - (AC - Circuits) 

Aufgabenstellung: 

Berechnung von passiven R-C-L-Wechselstromschaltungen unter Anwendung von komplexen Zah- 
len. 

Calculation of passive R-C-L-AC-Circuits using complex numbers. 

Die beschriebene Prefix-Notation laftt sich auch auf Berechnung von Wechselstromschaltungen uber- 
tragen. Als Beispiel dafur folgt die Vorstellung eines Listings zur Berechnung von gemischten 
R-C-L-Wechselstromkreisen, die sich auf Reihen- und Parallelschaltungen zuriickfuhren lassen. Da 
DERIVE die komplexen Zahlen, sowie die auf diesen Zahlenbereich erklarten Rechenoperationen fest 
implementiert hat, bietet sich bei der Berechnung charakteristischer GrofJen von Wechselstromkrei- 
sen die Anwendung komplexer Zahlen an. Im Rahmen dieses Aufsatzes beschrankt sich die Pro- 
grammierung auf die Berechnung des „komplexen Scheinwiderstands" 10 Z, dessen Betrag Zabs und 
des Phasenverschiebungswinkels a. 

Die Zeilen #1 bis #4 erklaren die (Grund-) Bestimmungsgleichungen fur den Scheinwiderstand einer 
Reihen- bzw. Parallelschaltung aus R und L bzw. C. Dabei ist ZE_ .... ein Element der komplexen 
Zahlen und wird in Normalform (a + T.b, mit f = - 1) dargestellt. Durch diese vier Grundgleichungen 
lassen sich alle weiteren, aus Reihen- bzw. Parallelschaltungen hervorgegangenen Wechselstrom- 
netzwerke symbolisch beschreiben und berechnen. Die Faktoren „10" 3 “ bzw. „10" 6 “ ermoglichen die 
Eingabe der Induktivitat in mH bzw. der Kapazitat in pF. 


1 0 

Unter diesem Terminus versteht man den ..Scheinwiderstand in symbolischer Darstellung unter Anwendung der komplexen 
Zahlen“. 







D 


-N-L#15 


H. Scheuermann & N. Wild: Elektrotechnik 2 


p 35 


-3 

#1: ZE_SER_RL(r , 1 , f) := r + i- 2 -tt- f ■ (1 -10 ) 

L 

ZE_SER_RC(r , c, f) := r 

# 2 : -6 

2-TT.f-Cc.10 ) 

1 


L 


-3 

2-tt- f- (1-10 ) 

1 

ZE_PAR_RC(r , c, f) := 

# 4 : 1 -6 

+ i.- (2-tt- f - (c-10 )) 

r 

Die folgenden Zeilen #5 und #6 legen die allgemeine Operationsvorschrift fest, wie der komplexe 
Scheinwiderstand einer Reihen- bzw. Parallelschaltung aus zwei komplexen Scheinwiderstanden zl 
und z2 gebildet werden kann. 

Expr. #5 und #6 provide how to build the complex impedance of a series- and parallel connection from two 
complex impedances zl and z2. 

#5: Z_5ER(zl, z2) := zl + z2 

zl- z2 

#6: Z_PAR(zl, z2) := 

zl + z2 

So erklart z.B. der Term z_SER(ZE_SER_RL(rl,l .f) ,ZE_PAR_RC(r2,c,f)) den komplexen Scheinwi- 
derstand der Schaltung bestehend aus R1 in Reihe zu L in Reihe zu R2 || C (vgl. Abb. 1 1 ) 

Zur Berechnung des Phasenverschiebungs- 
winkels a und des Betrages des 

Scheinwiderstandes werden die vorhandenen 
DERIVE-Funktionen PHASE bzw. ABS 

verwendet. Die Division durch das Symbol „° “ 
in Zeile #8 bewirkt die Ausgabe von a im 
GradmaG>. Der komplexe Gesamt- 
scheinwiderstand wird mit Z(x) bezeichnet. 

Dabei ist x als Platzhalter fur die moglichen 
Parameter eines Wechselstromkreises (zB bei 
der Darstellung von Ortskurven) erklart, also: 

X €E {R, C, L, f}. 



Abb 11: Schaltbild eines R-C-L-Wechselstrom- 
kreises 


ZE_PAR_RL(r, 1, f) := 

1 

#3: — 

r 


#7: 

Z(y') 

:= 

#8: 

a := - 

PHASE (ZOO) 

1° 



#9: 

zabs 

:= |ZC*)I 


Diese Programmierung wird zur Losung (vgl. Abb. 11) der nachfolgenden Aufgabe den weiteren 
DERIVE - Zeilen vorangestellt Oder im „Hintergrund“ aufgerufen. 
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Aufgabenstellung: 

Fur den in Abb. 12 dargestellten Wechselstromkreis sind die folgenden Werte gegeben: 


R 1 = 25Q; L = 79,5mH; 

R 2 = 12,5Q; C = 40jliF. 

a) Man bestimme den komplexen Scheinwi- 
derstand sowie dessen Betrag und Pha- 
senverschiebungswinkel a fur f = 50 Flz. 

b) Fur welche Frequenz f wird die Schaltung 
reell? 

c) Man zeichne die Widerstandsortskurve 
und parametrisiere diese fur 10Hz < f < 
1 00Hz mit einer Schrittweite von 1 0Hz 



Abb 12: Schaltbild laut Aufgabenstellung 


Durch die „Vorarbeiten (Zeilen #1 bis #9) gestaltet sich das Losungslisting (vgl. Abb. 13) zur vorgege- 
benen Aufgabe recht kurz. Zeile #12 stellt die symbolische Beschreibung der R-L-C- Wechselstrom- 
schaltung dar (mit x = f). Nach Festlegung der Frequenz f konnen durch Vereinfachung der Zeilen #8, 
#9 und #12 der komplexe Scheinwiderstand, dessen Betrag und Phasenverschiebungswinkel direkt 
abgerufen werden (vgl. Zeilen #14, #15 und #16). Die Zeile #19 erklart die Frequenz f (wieder) als 
Variable. Der Aufruf IM(Z(f)) bestimmt den Imaginarteil der Wechselstromschaltung, der zur Bestim- 
mung des reellen Scheinwiderstandes gleich Null gesetzt wird (Zeile #20, Losungen stellen die Zeilen 
#21 und #22 dar). Die Ortskurve wird uber die Berechnung der Wertepaare in Zeile #25 geplottet (der 
Realteil wird auf die I.Achse, der Imaginarteil auf die 2.Achse abgetragen). der VECTOR-Befehl in 
Zeile #26 erstellt die Wertetabelle in #27 und ermoglicht die Parametrisierung der Ortskurve laut Auf- 
gabenstellung. (#30: [ RE ( Z ( x ) ) , IM(Z(x))]) 
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#2S : Solution c) 

#26: [RE(Z(x)) , IM(Z(x))] 

#27: APPEND([[f , Re, Ini]], VECT0R([f, RE(Z(x)), IM(Z(x))], f, 10, 100, 10)) 


" f 

Re 

Im 

10 

25.4S54 

3.42153 

20 

27.0126 

6.81133 

30 

29.7965 

10.0526 

40 

34.1S39 

12.8077 

50 

40.5301 

14.2874 

60 

4S.72S1 

13.0062 

70 

57.1S15 

7.11733 

SO 

62.2729 

-3 .60468 

90 

61.0207 

-15.7024 

. 100 

54.6439 

-24.684 


|[RE(Z(x)), IM(Z(x))]| 


Abb 13b: Losung zur Aufgabe R-C-L-Wechselstromkreis 



Abb 13c: Darstellung der Ortskurve 


Problems to be solved: 

a) Find the complex impedance, its absolute value and the angle of phase displacement for f = 50Hz. 

b) Which frequence f will make the circuit real? 

c) Give the parametric form of the locus and plot it (valuetable in #28). 
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7 Einschaltvorgange - Turn - on - processes 
Aufgabenstellung: 

Gegeben ist die Parallelschaltung aus zwei Rei- 
henschaltungen R 1s C bzw. R 2 und L (vgl. Abb 
14). An den Klemmen liegt die Gleichspannung 
U. Gesucht ist der Zeitpunkt t fur den I.C = I.L 
gilt, mit 

R-i = 100 Qk, C = 1 nF 
R 2 = 20 Qk, L = 1 H. 


Abb 14: Einschaltvorgang 



# 1 : 

# 2 : 

#3: 

#4: 

#S: 

# 6 : 

#7: 

# 8 : 

#9: 

# 10 : 

# 11 : 

# 12 : 

#13: 


#14: 


#15: 


#16: 


# 17 : 


We rtzuwei's ungen - Assignment of values 

5 4-9 1 

rl := 10 , r2 := 2-10 , c := 10 , 1 := lj 

Besti'mmung der Zen tkonstanten t1 und t2 der Rei henschal tungen 

Determination of the series connections' constants of time t1 and t2 

1 

t 1 := rl-c, t2 := 

r2 

Funkti'onsglei'chungen fur die Einschaltvorgange 

Equation for the turn-on- pro cesses 

u - t/T2 

IL( t ) := (1 - e ) 

r2 

u - t/rl 

IC(t) := e 

rl 

Zur Berechnung des Zei'tpunktes fur Stnomglei'chhei't gilt: 

For calculating the moment of equality of current: 

IL(t) = IC(t) 

-5 -5 - 20000- t -5 - 10000- t 

B-10 -u-5-10 ■ u ■ e =10 ■ u ■ e 

- 20000- t - 10000- t 

1 e e 

20000 20000 100000 

soLve > Memory full!!! 

Substi tute t = -z/10M 

4 4 

- 20000- (- z/10 ) - 10000- (- z/10 ) 

1 e e 


20000 


20000 


100000 
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2-z 


#18: 


#19 : z = LN 


20000 20000 100000 
^ JlOl 1 ^ 

10 j 
1 


LN 


. 10 
JlOCL 
10 


10 


# 20 : 


# 21 : 


10 


IBM 


Dieses Beispiel zeigt, dass auch DERIVE (Version 2.52) kein Alleskonner ist. Wichtig erscheint uns, 
den Schulern auch solche Beispiele vorzustellen 11 . Dies verfolgt die Absicht, den Lernenden zu einer 
kritischen Haltung gegenuber dem Rechner bzw. der Software anzuregen und damit dem Entstehen 
einer bequemen und unreflektierten Rechnerglaubigkeit entgegenzuwirken. 

This example shows that DERIVE (v 2.52) isn't a jack-of-all-trades. For us it is important to give our students 
such examples, too. It is our aim to encourage a critical attitude towards the computer and the software. 

Later Derive versions have no problems in numerical solving the equation, but is not able to solve the 
equation in Exact Mode. 


N SOLVE 


1 


L 20000 


- 20000- 1 - 10000- 1 ^ 
e e. -6 

= , t = (t = 9.983407889-10 ) 

20000 100000 J 
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11 Diese Aufgabe aus dem Elektrotechnikunterricht wurde von den Schulern mit der Bitte an uns herangetragen, ihren geschei- 
terten Losungsversuch mit DERIVE zu untersuchen. Die Moglichkeit, daft das Programm bei dieser Aufgabe „aussteigt“ wurde 
nur sekundar in Erwagung gezogen. Primar erwarteten sie von uns einen Hinweis auf die „richtige“ Einstellung der Programm- 
parameter bzw. auf die Wahl der losungsbringenden Tastenkombination. 

(Students came to us with the wish to investigate their attempt to find a solution with DERIVE, because they had failed. It 
was interesting that they didn't consider DERIVE to fail. They expected us to give them some advice for the correct setting 
of any “options" or to show them the right choice of keystrokes.) 
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Titbits from Algebra and Number Theory (3) 

by Johann Wiesenbauer, Vienna 

One of the most important tools of number theory is the so-called Moebius-function (yes, the 
same Moebius who invented the famous strip!) which is defined on the set of natural numbers 
in the following way: 

. . f 0 if n is divisible by a square 

\i(n) := ( , 

{(-1) if n is a product of k distinct primes 

As for its implementation in DERIVE I have good news and bad news for you. First the bad 
ones: Unlike some other computer algebra systems there exists no built-in p- function in 
DERIVE, and what is more, it is rather cumbersome to write a program for it. Here is the re- 
sult of all my efforts to do so: 

LPD(n, s) := 

If NEXT_PRIME(n - 1) = n 
n 

ITERATE (I F(t A 2 > n, n, IF(M0D(n, t) = 0, t, NEXT_PRIME(t))) , t, s) 

(Returns the least grime divisor of n which is greater or equal to s > 1. In case, it doesn't ex- 
ist, the result will be n .) 

Examples: 

LPD(2431, 12) = 13 


FACTOR (2431) = 11-13-17 


f 

ITERATE 

f 

IF 

r 

x = 1, X, 



2 

v 




Xx 5 r x 

i . : 


x 3 

3 JJJ 


(Yields the list of all grime factors - not necessarily distinct - of n in increasing order.) 
Examples: 

PFLC2431) = [11, 13, 17] 

PFL(IOOO) = [2, 2, 2, 5, 5, 5] 

PREMOEBIUS(pl) :: 

If DIMENSION(pl) > 1 A FlCpIlk - pl|(k - 1), k, 2, DIMENSION(pl)) = 


(_l) A DIMENSION(pl) 


(One more auxiliary function, but . . .) 
p(n) := PREMOEBIUS(PFL(n)) 


(. . . here it is, at last! If somebody comes up with a better solution, I would be dead keen on 
hearing about it!) 






Johann Wiesenbauer: Titbits 3 


D 


-N-L#15 


p 41 


Examples: 

MC2431) = -1 
MC1000) = 0 
1X1122) = 1 

And now the good news: In the forthcoming version 3 of DERIVE all the code above can be 
replaced by one single line! 

M_(n) := n(IF(PRIME(IO, - 1 , 0), k, FACTORS ( FACTOR (n))) 


In this short program I have taken advantage of 4(!) new features of DERIVE3, namely: 

• The new fast routine PRIME (k) replaces the former inferior primality test 
NEXT_PRIME(k-l)=k. 

. The menu option FACTOR is now also available as a function. 

• The output of FACTOR(n) is transformed into a vector of the corresponding prime powers 
of n by FACTORS. (Unfortunately, the primes of n are still not directly accessible.) 

. The domain of the variable k in the product above may now be given as a vector. 

Examples: 

H_(2431) = 0 

p_(1000) = 0 
M_(1122) = 0 


Comparing with the results from above you will note that this routine does not work. So let's investi- 
gate what is wrong and how to adapt Johann's formula for recent versions! 

This was Derive 3 This is Derive 6: 


FACTOR (1000) 
3 3 


2 • 5 
FACTORS 
[ 1000 ] 
FACTORS 



1000 ) 

FACTOR (1000) ) 


3 3 

FACTOR CL00CO = 2-5 


FACTORS (1000) = 


2 3 
5 3 


So I “reinvent” an old_factors-function for 
DERIVE 6 


olcLfactors(n) VECTOR 
olcLfac tors (1000) = [8, 125] 


v , v, FACTORS(n) 
1 


p_(n) := ri(IF(PRIME(k) , -1, 0), k, old_factors(n)) 


M-C2431) = -1 

M_(1000) = 0 
M_(1122) = 1 
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This works, but an interesting happens simplifying this function: 

p_(n) := F[(IF(PRIME(k) , -1, 0), k, ol cLf actor s(n)) 

M-tn) :z 

If PRIME?(n) 

-1 

0 


[M-C2431), M-C1000), M-CH22)] = [0, 0, 0] 

Finally we will not keep secret, that the DERIVE versions from today provide a built-in 
Moebius-function MOEBlUS_MU(n): 

[MOEBIUS_MU(2431) , MOEBIUS_MU(1000) , H0EBIUS_MUC1122)] = [-1, 0, 1] 
Let's come back to Johann's Titbits from 1994. 


One of the most important applications of the p- function is the so-called Moebius inversion 
formula (cf. [1]). Without going into details I'd like to explain its use in connection with the 
factorization of polynomials of the form x n - 1 . If we define CJx) for any natural number d by 

c^n^-ip 

S|rf v ' 

(8 runs through all positive divisors of d or, what amounts to the same, through all squarefree 
positive divisors of d) then it turns out that this rational expression in x can always be simpli- 
fied to a polynomial with integer coefficients which is called the ct h cyclotomic polynomial. 
Moreover, the representation of x n - 1 by a product of irreducible polynomials over the field 
of rational numbers is given by 

x”-\ = Y\CAx). 

d\n 

And here is the DERIVE-implementation: 

ENLARGE (v, a) := 

If M0D(viDIM(\0, a) = 0 
v 

APPENDCv, a-v) 

(An auxiliary outer operation that is used for a recursive adding of new squarefree 
divisors to v) 

PRESFDL(pl) :z 
If DIM (pi) < 2 
APPEND([1] , pi) 

ENLARGE(PRESFDL(VECT0R(p1|k , k, DIM(pl) - 1)), p1|DIM(p1)) 

(First we need a list of all squarefree divisors of n ; this auxiliary function does the main job.) 
SFDL(n) :z PRESFDL(PFL(n)) 

SFDL(1000) = [1, 2, 5, 10] 

(That's it. Any idea, how a routine for all divisors of n could look like?) 
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About Cyclotomic Polynomials (JB) 

A cyclotomic polynomial of order d is a polynomial given by 

k = i 

with Q = the primitive d th Root of Unity in C. 

A number r is an n th Root of Unity if r n = 1 and it is a primitive n th Root of Unity if n is the smallest 
integer of k = 1, . . n for which / = 1. 
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Let's have an example: What is the cyclotomic polynomial of order 6 = <J> 6 (x)? 
We can use the function from page 9 to find all Roots of Unity of order 6. 

L-tt-2- k/n 


#1: 1 

NTH_R00T5_0 F_UNITY(rO := VECTOR (e 


i k. 

o, n - 1) 

#2: 

noot6 

:z 1 

NTH_R00T5_0 F_U N ITY (6 ) 




#3: 

nootS 


1 

1 . 

73- i 1 

J3.L 

_i 

1 -/3-i . 1 -/3-i . 


-L , "h 

2 

2 1 2 

2 

l -L i 

2 2 ' 2 2 

k 

#4: VECTOR (VECTOR (r , 

r , roots) , k , S) 





“I 

1 

J3-L 

1 73. i. 


1 

J3-L 1 73-d. 


_L 

2 

+ 

2 

2 2 


2 

2 2 2 


-1 


L 73- i. 

1 73. i. 


1 

73-d. 1 73-d. 


-L 


+ 

2 2 

2 2 

_L 

2 

2 2 2 

#S : 

1 


-1 

L J3- L 

i 

1 73. i. 

-1 

1 

1 -1 
43- i 1 J3- i 


-L 


2 2 

2 2 

_L 

2 

2 2 2 


1 

1 

J3.l 

1 73- i. 

-1 

1 

73-d. 1 73-d . 


_L 

2 

2 

2 2 

-L 

2 

"h 1 

2 2 2 


. i 


1 

1 

1 


1 1 

#S: 

(x - 

roots ) ■ (x - 

roots ]i = jS 

■ 




2 S 


Only root #2 and root#6 have k = n = 6 as smallest integer to give r k = r n = 1 . A table of cyclotomic 
polynomials shows ® 6 (x) = x 2 -x + l and CP(x , 6) from above gives the same result, of course. 

Try to find ® 4 (x) and ® 5 (x) in this way. What about ® 9 (x)? 

By the way, Johann mentioned the cyclotomic polynomial with d = 105, the next ones with a d of in- 
terest are with d= 385, 1365, 1785. Check it! 

Reference: CRC Concise Encyclopedia of Mathematics 


Paolo Ruffini 

Paolo Ruffini was an Italian mathematician and physician. He was born 1765 in 
Valentano and he died in 1822 in Modena. 

Paolo Ruffini was professor for Applied Mathematics and for Practical Medicin at the 
University of Modena. As a mathematician he contributed for the solution of algebraic 
equations. In 1799 he published his conjecture that general equations of order 5 can- 
not be solved algebraically but his proof was uncomplete. 
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Titbits(3) - then and now 

Johann Wiesenbauer, Vienna University of Technology, June 2007 

This is a revised version of my Titbits(3) in the DNL #15. Note that if there are two functions 
with the same name, the first one is usually the old version (sometimes with small notational 
changes), which is then overwritten by an updated version that makes full use of all the features 
of Of W 6.10 along with other improvements. 

The first function there, namely 

lpd(n, s) := 

If NEXT_PRIME(n - 1) = n 
#1: n 

ITERATE(IF(t A 2 > n, n, IF(M0D(n, t) = 0, t, NEXT_PRIME(t))) , t, s) 

was supposed to return the least prime divisor p of n that is not below a given integer bound s >1, 
but actually the output p is only an ordinary divisor of n (not necessarily prime!) with the prop- 
erty p>=s. In the following this is corrected, although this error had no implications as to the 
calls in the subsequent functions. After introducing additionally the default case s=2, this func- 
tion could nowadays look like this: 

lpd(n, s := 2, t_) := 

Prog 

If PRIME?(n) 

RETURN IF(n > s, n) 
t_ := s 

#2 : Loop 

If t_ A 2 > n 
RETURN ? 

If M0D(n , t_) = 0 
RETURN 1 pd(t_, s) 
t_ :+ 1 

100 

#3: lpd(10 + 9) = 3221 

The "dual" function to lpd(n,s) would be minprime(n,s) (by the way, an only slightly adapted ver- 
sion of the min_prime(n,s) in my recent update of the Titbits(2)) that computes the least prime 
divisor p <=s of n, if it exists. For the default cases of s both functions lpd(n) and minprime(n) 
should coincide, returning the least prime divisor of n, if it exists. 

mi nprime(n , s := a>, t_ := 2, d_ := 4) := 

Prog 

If PRIME?(n) 

RETURN IF(n < s, n) 

If s = oo 

s := FLOOR(Jn) 

Loop 

#4: If t_> s 

RETURN ? 

If M0D(n , t_) = 0 
RETURN t_ 

If t_ < 5 
t_ :+ 1 
t_ : + d_ 
d_ := 6 - d_ 
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#5: mi nprime(10 + 9) = 3221 

The next function pf l(n) was supposed to give a list of all (not necessarily distinct) prime factors 
of n in increasing order. 

( ( f x ( x 

2 2 

pfl(n) := ITERATE IF x = 1, x, APPEND(x , Tx 1), , lpd , 

#6: 2 1 L 3j x x 

l l L 3 v 3 


I, x, [[], n, 1 pd(n , 2)] 


#7: pfl (360) = [2, 2, 2, 3, 3, 5] 

Well, it works, as the example above shows, but meanwhile the much more powerful built-in func- 
tion factors(n) has become available. 


FACTORS (360) =32 


What I was actually aiming at in my article then was an implementation of the Moebius p- 
function. By definition p(n)=0, if n is not squarefree, i.e. at least one of its prime divisors has a 
multiplicity greater than 1, as the prime divisors 2 and 3 for n=360. On the other hand, if n is a 
product of r distinct prime numbers, then p(n) = (-l)'Y. Hence, my implementation of p(n) looked 
like this then 

premoebius(pl) := 

If DIM(pl) > 1 a n(pl|k - pi j. C k - 1), k, 2, DIM(pl)) = 0 
#9: 0 

(-1) A DIM(pl ) 

#10: p(n) := premoebi us(pfl (n)) 

#11: [pCD, H(30), p(60) ] = [1, -1, 0] 

Oh well, believe it or not, but even this clumsy implementation was a significant progress at the 
ancient times of Derive 2.xx! Needless to say that since the availibility of the crucial function 
factors(n) in version 3 of Derive this has been greatly simplified and now all this has become 
the following onelinerl! 

MOEBIUS_MU(n) := Il(- MAX(2 - v_ , 0) , v_, FACTORS(n)) 

#12: 2 


There are a number of very important applications of Moebius p-function, and as an example I 
gave an implementation of the so-called cyclotomic polynomials. Basically the n-th cyclotomic 
polynomial can be defined as the product of all linear polynomials x - exp(2-m k/n), where k runs 
through all numbers in l,2,..,n, which are coprime to n. 


Instead of listing all the clumsy routines for cyclotomic polynomials in the original Titbits(3), 
which are totally outdated by now anyway, let’s use Derive to compute §24(x) according to this 
definition 
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48 47 46 43 42 41 40 39 

#22: CYCLOTOMIC(105 , x) = x +x +x -x -x -2-x -x -x + 

36 35 34 33 32 31 28 26 24 22 20 17 

X + X + X + X + X +X -X -X -X -X -X + X 

16 15 14 13 12 9 8 7 6 5 2 

+ X + X + X + X + X -X-X - 2-X - X - X + X + X + 1 

80 79 78 75 74 73 69 68 67 

#23: CYCL0T0MIC(165 , x) = x +x +x -x -x -x -x -x -x 

65 64 63 62 60 59 58 54 53 52 50 

+ x + 2-x + 2-x + x -x -x -x -x -x -x +x 

49 48 47 46 44 43 42 41 40 39 

+ 2-x + 2-x + 2-x + x -x -x -x -x -x -x 

38 37 36 34 33 32 31 30 28 27 26 

x -x -x + x + 2-x + 2-x + 2-x +x -x -x -x 

22 21 20 18 17 16 15 13 12 11 7 

-x -x -x + x + 2-x + 2-x +x -x -x -x - x - 

6 5 2 

X - X + X + X + 1 

96 95 94 91 90 89 83 82 80 

#24: CYCL0T0MIC(195, x) = x +x +x -x -x -x -x -x +x 

79 78 77 75 74 68 67 65 64 63 62 
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6 5 2 

X - X + X + X + 1 

No, of course, this is not yet the solution, rather sort of a red herring. 
But here it comes at last! 

#25: SELECT(|n(SUBST(TERMS(CYCLOTOMIC(n, x)), x, 1))| > 1, n, 200) 


#26: 


[105, 165, 195] 





DERIVE POETRY 


p 50 


D-N-L#15 


The spirit of the Plymouth DERIVE conference inspired two of our DUG-members to write 
poems about DERIVE and its symbol - the little gecko. 

Sergey Biryukov is not only a skilled DERIVE user but also a poet. Like a great actor he re- 
cited his poem and for all of us who were present this was a great moment: 


DERIVE - ing physics at the class 
We use the air and the glass 
We Plot the rays, we mirror Plot 
Light source we represent as Dot. 

We use Expand and Simplify 
We say hand calculus good-bye 
DERIVE we use as tool for mind 
Is a unique software kind. 

It gives numeric, graphic face 
Its vectors handle every space 
An analytic face - the main 
It is intelligence retrain. 

If use DERIVE in mathes need 
You will undoubtedly succeed! 


Jim Hall wrote a little song, too: 


I'M A LITTLE GECKO 

(to the tune of 'I'm a Little Tea Pot') 

I'm a little gecko 

Slim, not stout. 

Here is my tail and 

Here is my snout. 

If you pull my tail off, 

Hear me shout: 

Ha - ha - ha - it will 

Grow back out!' 


